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Abstract 

Transitions from classical to quantum behaviour in a spin system with two 
degenerate ground states separated by twin energy barriers which are asymmetric 
due to an applied magnetic field are investigated. It is shown that these tran- 
sitions can be interpreted as first- or second-order phase transitions depending 
on the anisotropy and magnetic parameters defining the system in an effective 
Lagrangian description. 
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1 Introduction 



Barrier penetration by tunneling processes is a purely quantum phenomenon 
which does not arise in classical physics where only processes leading over the 
barrier, e.g. thermal activity according to a classical Boltzmann distribution, 
yield a nonzero barrier transition rate. At finite temperature, either tunneling 
from thermally excited states ( "temperature assisted tunneling") or thermal fiuc- 
tuations over the barrier ("thermal activity") dominate the transition rate, and 
the crossover from temperature assisted tunneling to thermal activity can be un- 
derstood as a phase transition from the quantum phase to the classical phase of 
a physical system which is of either first or second order. 

Whereas the general theory of these phase transitions in an abstract potential 
barrier setting is well known and clearly understood (there is a remarkable 
similarity to the Maxwell theory of phase transitions in the Van der Waals gas), 
only very few models are known which allow an explicit and analytic investiga- 
tion of the phase transitions in decay and transition rates which may even be 
accessible to experimental verification. Hence the recent discovery that spin sys- 
tems provide examples which exhibit first- and second-order phase transitions 
aroused interest in the investigation of such systems. In particular, a large 
spin in an XOY easy plane anisotropy with easy y-axis can be shown to exhibit 
both first- and second-order phase transitions depending on the value of the 
anisotropy parameter. 

In the following, we consider this spin system with an additional applied mag- 
netic field and investigate its infiuence on the dominant transition process. We 
begin with the presentation of the model and its effective semiclassical Lagrangian 
in Section |^, and then review the theory of temperature assisted quantum tunnel- 
ing and thermal activity in Section ^. Section ^ contains some analytical results 
which guided the numerical calculations presented in Section ^ and discussed in 
the concluding Section ^. 

2 The model and its semiclassical approxima- 
tion 

We consider a giant spin in an XOY easy plane anisotropy with easy axis along 
the x-direction and external magnetic field B in the y-direction, perpendicular to 
the easy direction. The corresponding Hamiltonian H involves the spin operator 

S and is given by 0, || 

H = K{Sl + \Sl) - 2nBBSy (1) 

where the easy xy-plane demands A < 1. 

The physical situation described by this Hamiltonian is illustrated in Fig. |l] 
where the spin operator is represented as a classical spin vector G S^. For zero 
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Figure 1: Classical visualization of the XOy-easy plane x-easy axis spin system 
with applied magnetic field in the y-direction. 

magnetic field B = the ground state is twofold degenerate, the classical spin 
vector pointing along the positive or negative x-direction, i.e. along the easy axis. 

Under the infiuence of an applied magnetic field in the ?/-direction, there are 
still two degenerate spin ground state directions S'^]„, S'^]„ in the easy x?/-plane 
moving towards the y-direction with increasing field B. 

To change the direction of the spin from one of these ground state directions 
to a neighbouring one, one has to avercome an energy barrier, moving the spin 
along either path S or path L. 

To study quantum tunneling and classical thermal effects of the discrete spin 
system described by H, we convert the spin operators to a continuous potential 
problem. This can be achieved with the help of spin coherent state path integrals 
or using the Villain transformation . Both approaches yield a semiclassical 
description of the quantum system given by the effective Lagrangian 

/:(0,0) = lM(0)02-y(0) (2) 

where 

V{(f)) = KXs'^{sin^(l)-a)'^ (3) 

^^'^^"2K(l-Asin2 + aAsin(/))' " KA^' 
Here may be interpreted as a spherical parameter of the classical spin vector 

s = s(sin 6* cos 0, sin 6* sin 0, cos 6*). (5) 
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The semiclassical approximation is exact in the hmit of large spin, s — oo, in 
the entire range of the anisotropy parameter A, < A < 1 . 

The shape of the potential V{(p) is shown in Fig. ^, together with the field- 
dependent mass M{(p) which is a special feature of this model. For small magnetic 
and anisotropy parameters, i.e. A <^ 1, a <^ 1, one can approximate the mass 
function by a constant value M[(f)) ^ ^ [|], but we are here particularly in- 
terested in the effect of nonconstant mass on quantum tunneling and thermal 
activity. Therefore we restrict the mass only to be positive which yields the con- 
dition A(l + a) < 1 for the magnetic parameter a, but keep the full ^-dependence 
of M(0). 

The degenerate spin ground states are given by the two different types of 
minima of V{<f)) at 2ln + arcsina and {21 + 1)tt — arcsina. These minima are 
separated by a small barrier S with height 

v(^^'^=KXs\l-ar=:Es (6) 

and a large barrier L with height 

v(^^^=KXs{l + af=:EL (7) 

which correspond to the paths S and L shown in Fig. |1|. 

Since and (j) + 2% describe the same physical state, we restrict ourselves to 
the first twin barrier pair at (^5 = | and (fi = The maximum positions (^s 
and (fL are called "sphalerons" in the usual field theoretical terminology. The 
vacua surrounding these barriers are denoted by 

$5 ■* = arcsin a ~ <i>^'' = 27r + arcsin a (8) 

and 

$5 ^ = TT - arcsin a = . (9) 

We observe that barriers of different height only exist for < a < 1 . For vanishing 
magnetic field (a = 0), the two barriers are equally high ||^, whereas in the limit 
a — i> 1, the magnetic field dominates the easy axis effect and there is only one 
ground state pointing along the y-direction. 

3 The theory of quantum tunnehng and thermal 
activity 

We consider transitions between spin states built around the two degenerate vacua 
at finite temperature, i.e. we assume the quantum spin states to be populated 
according to a Boltzmann distribution. 
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The rate of thermal activity over the small or large barrier is to first order 
given by the Boltzmann factor, 

f,(r)~e^-(^), ~S^{T):=^ (10) 

with i = S,L and = 1- Ss{T) and Sl(T) are called the thermodynamic 
actions of the small and the large barrier, respectively. 

The temperature assisted tunneling rate can be estimated by a Boltzmann av- 
erage over the tunneling probabilities from excited states with energy E. These 
tunneling probabilities can be approximated by the semiclassical WKB expo- 
nents, Pi{E) = e-^^^^\ 

W,{E) = 2V2 /*; ^""^ ^M(0)(y(0)-i?)d0, z = S,L (11) 

where $5 (E), $2.' \E) are the turning points for the small (5") or large (L) 
barrier at energy E, i.e. the solutions of the equation V{(f)) = E. It is easy to see 
that 




and 



$«(E) = 

^l\e) = TT — arcsin | a 
= 27r + arcsin(a 



< ^ < ^5 (12) 



E 



KXs^ 



0<E<El. (13) 



E 



KXs^ 

Taking the Boltzmann average over the tunneling probabilities from excited 
states Pi{E) yields the temperature assisted tunneling rate, 

TiiT) = / dEe~^Pi{E) = / dEe-^-^^'^^\ (14) 
JO Jo 

This integral can be estimated by the steepest descent method, using the concept 

of periodic instantons [0, ^ . These are classical solutions of the Euclidean Euler- 

Lagrange equations of the semiclassical Lagrangian (^, i.e. 

(dots now denote derivatives with respect to Euclidean time r = it) with finite 
energy E as integration constant in the first integral, 

^M(0)02 - y(0) = -E. (16) 
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For < E < Es or < E < El there are solutions v^f (t), ^^ij) oscillating 
around the small or large barrier with period ps{E), pl{E), respectively. These 
solutions can be derived by integrating eq. (|1^) which leads to elliptic integrals, 
but it is impossible to solve the resulting expressions for the explicit r-dependence 
of the solutions ^pg, ipf which are visualized in Fig. ^ 

Nonetheless it is possible to compute the periods and the Euclidean actions 
of these periodic instantons from eq. ( |l6l) which yields 



and 



Pi 



S^ = 2V2 ' JM(0)(y(0) - E)d<P + p,E = W,{E)+p,E. (18) 

In the steepest descent approach, the integral ( p!^ is dominated by the con- 
figuration satisfying 

^=P^{E), (19) 

i.e. the period of the periodic instanton has to be identified with the inverse 
temperature. This yields the usual periodic instanton tree approximation for the 
temperature assisted tunnehng rate, 

Ti{T) = e-^'(^^ (20) 

where Si(T) is the Euclidean action of the periodic instanton with period j; = pi. 

Hence, there are two different physical processes and two different energy 
barriers involved in the evaluation of the finite temperature spin transition rate. 
Ignoring the effect of the field dependent mass, it is obvious that the small barrier 
processes always dominate large barrier ones, and it has to be checked whether 
this changes by taking into account the field dependence of the mass, depending 
on the parameters A and a. Can transitions involving the large barrier become 
dominant over those involving the small barrier for specific values of these pa- 
rameters? 

Another question to be analyzed is that of the crossover from temperature as- 
sisted tunneling to thermal activity for the the small and the large barriers which 
can be understood as phase transitions of either first or second kind, depending 
on the shape of the function Pi{E) [|l], [1^]. This crossover can be visualized in 
a diagram showing both the thermodynamic and the periodic instantons action 
depending on temperature, i.e. T i— > {Si{T), Si{T)}. The phase transition occurs 
where the two curves intersect (sharp crossover, first-order phase transition) or 
join (smooth crossover, second-order phase transition) at lowest action. 

From eq. (p^, we obtain E = ^ and thus Si = J Edpi. The period Pi{E) 
of the periodic instantons usually decreases monotonically for increasing E near 



V2 



\ V(0) - E 



M(0) 
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E = 0. Hence if Pi{E) increases again after a certain critical value with increasing 
energy E ( "first-order behaviour"), the inverse function E{pi) is double-valued 
and so is Si{T). This leads to an intersection of the lower branch of Si{T) with 
Si(T) which is the first -order phase transition at temperature Tf, whereas the 
upper branch of Si(T) joins Si{T) at some temperature T*^, T^^^ < Tf. 

li Pi{E) is monotonically decreasing ("second-order behaviour"), there is only 
one branch of Si{T) which smoothly joins Si{T) at a temperature Tf = T/^; this 
yields a second-order phase transition. 



4 Analytical results 

To investigate the two questions mentioned, we calculate the functions ps{E), 
Ss(T) and pl{E), Sl{E) numerically and analyse their dependence on A and a. 
To obtain some analytical hints for this numerical analysis, we first discuss the 
E ^ 0, E ^ Ei limits of the periodic instanton periods Pi{E), i = S, L. 

In the limit E —>■ 0, the periodic instantons reduce to the usual (vacuum) 
instantons describing ground state tunneling at zero temperature through the 
small or the large barrier It is a special feature of this model (and an effect of 
the field-dependent mass) that although the vacuum instantons are not periodic, 
they reach the vacua between which they interpolate at finite time, i.e. Pi{E = 
0) < oo. Nonetheless, since Pi{E = 0) is very large, the vacuum instantons 
dominate the integral (|l^ also for T = and hence describe vacuum tunneling 
as can also be seen from the comparison with the vacuum WKB tunneling rate. 

Both Pi{E = 0) and the Euclidean action of the vacuum instanton, 

= SME = 0)) = V2 ,jM{(f))V{<t))dcP, (21) 
can be estimated explicitly in terms of elliptic integrals. Defining the parameters 



. o o 8JA(4 + Aa2) 
uo^ = AK^\s\ e = ^— (22) 



where 

2 471 2 471 2 2(a„ - 2v/A) 2 2(a+ - 2^) 

= = ^"^^ = (l-a)a^ ' ^"^^ = (l + «)a_ 

(23) 

and 

a+ = V4 + Aa2 + 2\f\ + a VA, a„ = V4 + Aa^ + 2\f\ - a\f\, (24) 
the vacuum limits of the periods are given by 

Ps{E = Q) = — ^ ( - (ai)^j n f arcsin (—) , {a^f , k) 

a;o(a- - 2vA)yar I L ^ V xa^J J 
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(ai)^F ^arcsin ^ — ^ , | 



Pl{E = 0) 



2./a 



{a^Y — (02)^ n ( arcsin ( — ) , (04)^, k 



+ (02) F ( arcsin 
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whereas the vacuum instanton actions are given by 



4s 



Sr 



< n arcsm — , (as) , k 

+ — — ^ arcsin — , A; U 

^a+a_ V V ^4 / /J 

4s { , — n arcsin — , {a^\ k 

a+al V \aj J 

(a^ -2^/X)a+ + 2^/Xaa_^ f . fl\ ^ 
H , F arcsin — ] ,k 



(25) 



(26) 



(27) 



(28) 



From these results, it is easy to check that the large barrier vacuum instanton 
action is always greater than that of the small barrier, Sl > Ss for all values 
of A, a considered. If both the large and the small barrier have second-order 
type periods pl{E), ps{E), then both Sl(T) and Ss(T) have only one branch 
which is strictly decreasing with increasing temperature T. These branches end 
at T = 7;^'^ with SiiT^'^) = SiiT^^). Eq. (|10|) yields Sl{T^) > Ss{T^) for all 
values of A, a leading to second-order behaviour for both barriers, hence 



Sl{T) > Ss{T), 



(29) 



i.e. tunneling through the small barrier always dominates tunneling through the 
large barrier if the crossover to thermal activity is a second-order phase transition 
for both barriers. Whether a first-order transition behaviour for one or both of 
the barriers allows large barrier tunneling to become dominant over small barrier 
tunneling has to be analyzed numerically. 

In the second limit E —>■ Ei, i = S, L, the periodic instantons reduce to the 
sphalerons 



E^Es 71" 

e^El Stt 



(30) 
(31) 



Near the maximum energy, E — Ei <^ 1, the periodic instantons can be approx- 
imated by small oscillations near the bottom of the inverted potential, and the 



9 



frequencies Ui of these oscillations determine the periods pi{E = Ei) = ^ of the 
static limits of the periodic instantons. To estimate these frequencies, one inserts 



(pi + Scj) into the Euler-Lagrange equation ([T5| ) and expands to second order 



in 5(f). This yields harmonic oscillator equations with frequencies 



ojs = ^j4:KXs^{l - a){l - X{1 - a)) (32) 



cj^ = ^AKXs^{l + a){l- X{l + a)). (33) 

Hence, the periodic instanton action curves Si{T) smoothly join the thermody- 
namic action Si{T) at Tf'^ = p^^E=E,) = t^- noting that T^^ > 
for A < |, but Tl^ < Tg^ for A > |. This suggests to investigate the parameter 
ranges AG (o, |) and A G (|, l) separately. 



5 Numerical results 

For vanishing magnetic field, a = 0, both barriers are equally high and the peri- 
odicity functions coincide, Pl{E) = ps{E) =: p{E), a situation already discussed 
[||. For < A < |, p{E) has second-order behaviour, whereas for | < A < 1, 
p{E) changes to first-order behaviour. 

A* = I remains a critical value of the anisotropy parameter if an applied 
magnetic field is considered. The main influence of the magnetic field on the type 
of phase transition is shown in Fig. |^ for = 1000, K = 1. For A = 0.45 < A*, 
Pl{E) changes from second-order behaviour to first-order behaviour when a is 
increased (Fig. ^L)), whereas Ps{E) which is not plotted is of second-order type 
for all values a G (0, 1). On the other hand, for A = 0.8 > A*, Ps{E) varies 
from first-order behaviour to second-order behaviour with increasing a (Fig. 
|^(S)), and one should note that the allowed values of a for A > ^ are restricted 
by A(l -|- a) < 1. Pl{E) for A = 0.8 which is not plotted exhibits first-order 
behaviour for all allowed values of a. 

These particular examples A = 0.45 and A = 0.8 are typical for the a- 
dependence of the type of transitions in the anisotropy parameter ranges A G 
(0, A*) and A G (A*,l). We can thus distinguish four different situations in the 
phase transition behaviour of the asymmetric twin barrier problem which are 
shown in Figs. § to ^ for values of A and a which yield clear shapes of the func- 
tions considered, again using = 1000, K = 1. These possible types of transition 
process combinations are summarized in Table |I]. 

For A G (0, A*), we have second-order phase transitions at both barriers (Fig. 
1^) for a < a* (A), or second-order phase transitions at the small and and first- 
order phase transitions at the large barrier (Fig. ^) for a > a* (A). The function 
a*{X) can be estimated numerically from the (A, Q;)-dependence of the period 
Pl{E). In Table @(L), some typical values of this critical parameter are given, 
showing that a* (A) decreases with A. For A < 0.25, there is no critical value 
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large barrier 



A = 0.45 





Figure 3: The periods of the periodic instantons of the large barrier pi for A = 
0.45 and of the small barrier ps for A = 0.8, plotted against e = for several 
values of a with = 1000, K ^ 1. 



11 



small barrier 



A = 0.3 a = 0.5 



large barrier 





Figure 4: Phase transitions for A = 0.3, a = 0.5 with = 1000, K — 1. (S-a), 
(L-a): The periods of the periodic instantons of the small and the large barrier 
plotted against e = ;^f;2- (S-a), (L-a): Periodic instanton actions Ss{T), Sl{T) 
and thermodynamical actions Ss{T), Sl{T) for the small and the large barrier. 
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small barrier 



A = 0.45 a = 0.8 



large barrier 



PS OA 
0.7p 
0.6 ■ 
0.5 - 
0.4 ■ 
0.3 - 



0.005 0.01 0.015 0.02 0.025 0.03 0.035 



(S-a) 



if As2 




(L-a) 



if As2 



S 9 




(S-b) 



T 




(L-b) 



4.5 



T 



Figure 5: Phase transitions for A = 0.45, a — 0.8 with — 1000, K — 1. 

(S-a), (L-a): The periods of the periodic instantons of the small and the large 
barrier plotted against e = -j^;^- (S-a), (L-a): Periodic instanton actions Ss{T), 
Sl{T) and thermodynamical actions Ss{T), Sl{T) for the small and the large 
barrier. 
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small barrier 



A = 0.7 a = 0.3 



large barrier 





Figure 6: Phase transitions for A = 0.7, a — 0.3 with — 1000, K — 1. 

(S-a), (L-a): The periods of the periodic instantons of the small and the large 
barrier plotted against e = -j^;^- (S-a), (L-a): Periodic instanton actions Ss{T), 
Sl{T) and thermodynamical actions Ss{T), Sl{T) for the small and the large 
barrier. 
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small barrier 



A = 0.9 a = 0.1 



large barrier 



PS 0.3 



0.28 



1 0.22 



0.2 



0.4 



(S-a) 



0.6 



if As2 



PL 



0.5 
0.55 

0.5 
0.45 

0.4 
. 35 



1 



0.2 0.4 



0.6 O.S 



(L-a) 



1.2 

E 



KXs^ 




s 




(S-b) 



(L-b) 



Figure 7: Phase transitions for A = 0.9, a — 0.1 with — 1000, K — 1. 

(S-a), (L-a): The periods of the periodic instantons of the small and the large 
barrier plotted against e = j^^- (S-a), (L-a): Periodic instanton actions Ss{T), 
Sl{T) and thermodynamical actions Ss{T), Sl{T) for the small and the large 
barrier. 
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A < 
small barrier 


• A* 
large barrier 


A > 
small barrier 


A* 

large barrier 


a < a*{\) 


second-order 


second-order 


first-order 


first-order 


a > a*(A) 


second-order 


first-order 


second-order 


first-order 



Table 1: Phase transition type combinations for the asymmtric twin barrier with 
field-dependent mass. 



Large Barrier 




Small Barrier 


A 


a*(A) 




A 


a*(A) 


< 0.25 






0.65 


0.06 


0.30 


0.85 




0.70 


0.12 


0.35 


0.63 




0.75 


0.17 


0.40 


0.45 




0.80 


0.23 


0.45 


0.31 




> 0.85 





(L) (S) 

Table 2: (L) Critical values a* (A) of the magnetic parameter at the large barrier 
for A < A*, (S) critical values a*(A) of the magnetic parameter at the small barrier 
for A > A* 

a*{X) < 1, i.e. the phase transition at the large barrier is of second order regardless 
of the magnetic parameter if the anisotropy parameter is sufficiently small. 

The second range of the anisotropy parameter, A G (A*,l), allows second- 
order transitions at the small and first-order transitions at the large barrier (Fig. 
1^) for a < a* (A), or first-order transitions at both barriers (Fig. |^) for a > 
a* (A). Some typical values of a* (A), now estimated numerically from the (A, a)- 
dependence of the period ps{E), are shown in Table 0(S). Here the critical value 
of the magnetic parameter increases with increasing anisotropy parameter. For 
A > 0.85, there is no critical value of a in the region (^0, ^ ~ 1 j restricted by the 
requirement of positive mass. For A < 0.60, the numerical calculations failed due 
to problems with the end-point integrations. 

We note that it is not possible to have first-order transitions at the small 
and second-order transitions at the large barrier for any allowed values of the 
parameters A, a. 

Moreover, the numerical analysis shows that processes involving the small 
barrier always dominate over those involving the large barrier even if one or both 
of the barriers exhibit a first-order phase transition behaviour. 
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6 Summary and conclusions 



Above we have analyzed the crossover from temperature assisted tunnehng to 
thermal activity for asymmetric twin barriers in a model with field-dependent 
mass describing a large spin in an XOF-easy plane with x-easy axis anisotropy 
and an applied magnetic field in the y-direction. 

The corresponding analytical and numerical analysis was guided by two ques- 
tions: 

1. Does the field-dependence of the mass allow the tunneling and/or ther- 
mal processes involving the large barrier to become dominant over those 
involving the small barrier? 

and 

2. What types of phase transitions are the crossovers from temperature as- 
sisted tunneling to thermal activity, depending on the anisotropy parameter 
A and the magnetic parameter a? 

Summarizing, the first question must be denied, i.e. small barrier processes 
always dominate large barrier processes which is physically obvious for constant 
mass from the shape of the potential and remains true even if the field-dependence 
of the mass is taken into account. 

But concerning the second question, the field dependence of the mass is of 
great importance because it leads to first-order phase transitons besides the usual 
second-order behaviour. This can already be observed for zero magnetic field ||^. 
Three of the four possible type of combinations of phase transition types were 
found: First-order transitions at both barriers, second-order transitions at both 
barriers and second-order transition at the small, first-order transition at the 
large barrier. The fourth possibility, first-order transitions at the small and 
second-order transition at the large barrier, did not arise. 

The types of combination of phase transitions depend on the values of the 
parameters A, a where A* = |, the critical value at which second-order behaviour 
turns to first-order behaviour for vanishing magnetic field (i.e. with equally high 
barriers), remains a critical value if a magnetic field is applied. For A < A*, the 
small barrier exhibits only second-order phase transitions; for A > A*, the large 
barrier exhibits only first-order phase transitions. The transition order at the 
other barrier, respectively, depends on the value of the magnetic parameter a 
and changes from second-order to first order at the large barrier at a* (A) for 
A < A* and from first-order to second-order at the small barrier at a* (A) for 
A > A* for increasing a, < a < min |l, — l|. 

We considered only the tree approximation of the tunneling rate to analyze 
its crossover to thermal activity. Taking one-loop corrections into account, i.e. 
calculating the fluctuation determinant prefactor in eq. (|1^), might perhaps 
smoothen a sharp intersection of the two curves. 
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Nonetheless, in experimental results a crossover between first and second order 
transitions can be observed, e.g., in molecular nanomagnets of spin 10-20, hence 
higher-order corrections are not expected to change the crossover behaviour sig- 
nificantly. The results derived here for the two-anisotropy model which is of high 
generality in small particle magnetism should therefore be helpful in experimental 
tests. 

Acknowledgements 

D.K. Park acknowledges support of the Deutsche Forschungsgemeinschaft (DFG). 

References 

[1] E. M. Chudnovsky, Phys. Rev. A46(1992)8011. 

[2] E.M. Chudnovsky and D. A. Garanin, Phys. Rev. Lett. 79(1997)4469. 

[3] M. Enz and R. Schilling, J. Phys. C19(1986)1765; C19(1986)L711. 

[4] J.-Q. Liang, H.J.W. Miiller-Kirsten, Jian-Ge Zhou, F. Zimmerschied and 
F.-C. Pu, Phys. Lett. B393(1997)368; J.~Q. Liang, H.J.W. Miiller-Kirsten, 
A.V Shurgaia and F. Zimmerschied, Phys. Lett. A237(1998)169. 

[5] E.H. Lieb, Commun. Math. Phys. 31(1973)327. 

[6] J. Villain, J. Phys. 35(1974)27. 

[7] E.M. Chudnovsky and L. Gunther, Phys. Rev. Lett. 60(1988)661. 

[8] J.-Q. Liang, H.J.W. Miiller-Kirsten, D.K. Park and F. Zimmerschied, Peri- 
odic Instantons and Phase Transitions in Spin Systems, Preprint. 

[9] J.-Q. Liang and H. J. W. Miiller-Kirsten, Phys. Rev. D51(1995)718; Phys. 
Rev. D46(1992)4685.. 

[10] A.N. Kuznetsov and P.G. Tinyakov, Phys. Lett. B406(1997)76. 

[11] H. Grabert and U. Weiss, Phys. Rev. Lett. 53(1984)1787. 



18 



